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We present a numerical computation of matrix elements of AI = 3/2 K tttv decays by using Wilson fermions. 
In order to extrapolate to the physical point we work at unphysical kinematics and we resort to Chiral Perturbation 
Theory at the next-to-leading order. In particular we explain the case of the electroweak penguins Ot.s which can 
contribute significantly in the theoretical prediction of e'/e. The study is done at P — 6.0 on a 24^ x 64 lattice. 



1. INTRODUCTION 

The study of non-leptonic kaon decays is a chal- 
lenging but important task, which is underlined 
by the recent accurate measurement of a non-zero 
e'/e parameter [1] (the first confirmed observation 
of direct CP violation) and the long-standing puz- 
zle of the A/ = 1/2 rule. 

For example, the electroweak penguins (EWPs) 
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although suppressed by the factor acm/ots when 
compared to the QCD penguins, are enhanced by 
the AJ =1/2 rule and tends to cancel signifi- 
cantly the contribution of the QCD penguins in 
the theoretical prediction of e'/e. 

Here we present the first lattice determination 
of the A/ = 3/2 i^T ^ tttt matrix elements (MEs) 
of the EWPs obtained directly with a two-pion 
final state. The strategy, explained in [2,3], con- 
sists in computing A^[I„**' =1=2 {t^''^\0-! ,s,\K^) at 
some unphysical kinematics on the lattice (called 
SPQR kinematics) and then using next-to-leading 
order (NLO) Chiral Perturbation Theory (xPT) 
to obtain the MEs at the physical point. 



2. OPERATOR MATCHING 

We now discuss the problem of the matching 
of lattice-regularised operators to the continuum 
renormalization scheme in which the Wilson co- 
efficients are calculated. Even if chirality is pre- 
served by the regularization, Oy^^ and Og^^ (the 
A/ = 3/2 parts of Oy^s) mix between themselves. 
In a Wilson-like regularization there is in general 
also mixing with operators of different naive chi- 
ralities. In the A/ = 3/2 case, SU{2) isospin 
symmetry forbids mixing with lower dimensional 
operators while CVS symmetry implies that mix- 
ing with other dimension six operators with dif- 
ferent chiral properties occurs only in the parity 
conserving sector [5]. The mixing is thus: 



Z{fj,a) 



(1) 
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where /i is the renormalization scale. With an 
obvious notation we define 2^77 (/ia), ZjsifJ-o.), 
2s7{fJ-a.), 2ss{fJ-a) to be the MEs of Z{fia). 

To compute Z{^a) we use the (mass inde- 
pendent) RI-MOM renormalization scheme which 
can be easily implemented in a non-perturbative 
way, also in the case of four fermion opera- 
tors [5]. Spurious effects due to the Goldstone 
boson contamination have been removed by us- 
ing the method explained in [6,2]. 

Once computed Z(fj,a) non-perturbatively, we 
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Figure 1. 2^^, , from the non-perturbatively determined Z{ii) as defined in Eq. 2. 

Table 1 

Results in GeV'' for the MEs at the physical 
point (the operators are renormalized in MS-NDR 
at /i = 2 GeV) obtained by using only points with 



have to check to what extent it follows the renor- 

malization group behaviour predicted by pertur- 
bation theory at the NLO. We thus compute 



with 



(2) 



where J has been computed in [7]. Z^*^^ should 

not depend on the scale /i (see Fig. 1). 

The observation that no clear 0(a'^p'^) lattice 
artifacts are visible in Zn, Z^j together with the 
very large anomalous dimension of Z^g,, Z^s (6 
and 16 at LO) suggest that higher orders in PT 
could contribute significantly in particular to the 
RG evolution of Zss- 

The central value for the renormalization con- 
stants at a given scale n is extracted in the fol- 
lowing way: we compute Z^^^ by fitting the 
plateaux to a constant in the interval 2.2 GeV < 

< 2.6 GeV. We then run Z^*^^ down to the de- 
sired scale /i. The systematic error is estimated 
by computing the variation of the value obtained 
at this scale from the evolution of the original val- 
ues (those used to build the combination of Eq. 2 
and Fig. 1) at the different scales /io- 

3. USE OF THE CHIRAL EXPANSION 

The NLO contribution in xPT, both in the 
quenched approximation and in the full-QCD 
case, has been computed for the SPQR kinemat- 
ics in [4]. One can thus think of extracting the 
low energy constants (LECs) of the LO and NLO 
(in the case of the EWPs 0{p^) and 0(j;2) respec- 
tively) from lattice data and to use them in the 
NLO formula for the physical kinematics. Both 
forms can be derived from the formula in the most 



Mk,Et^ < A. In the quenched case the values 
a = 0.1 and nio = 0.5 GeV have been used. 



fit A 




/Oc,)'"^^'' X 

{'^8)1=2 H „ f 


poly. 1.2 
poly. 0.8 


0.120(11) 0.023 
0.125(13) 0.012 


0.714(62) 0.031 
0.715(75) 0.037 


xPT 1.2 
xPT 0.8 


0.143(13) 0.817 
0.123(13) 0.482 


0.836(73) 0.815 
0.700(72) 0.626 


qxPT 1.2 
qxPT 0.8 


0.254(25) 13.2 
0.212(24) 5.95 


1.49(16) 10.8 
1.19(14) 4.91 



general kinematics (energy momentum injection 
via the weak operator allowed but with all the 
particles on-shell) which reads 

(MZfr^^ = A^{1 + [x logs],,„) + Mi + 
B^M^ + B^{p„+ + p,o) • pK + Bi{p^+ ■ p^o) (3) 

where we have symmetrized over Pt^+jPtt" since 
the two pion state is a pure 1 = 2 state. 

The SPQR kinematics corresponds to put the 
K and one tt at rest while the other tt is either 
at rest or has the minimum momentum allowed 
on the lattice. The previous formula can thus be 
re-expressed in terms of Mk, and S^, where 
i?7r is the energy of the (possibly) moving pion. 

Besides the fit with the complete form pre- 
dicted by xPT (Eq- 3), we also try a fit with the 
polynomial part only (the coefficients of which are 
now "effective" LECs) which we caU {M'J,^'> . 
Tab. 1 reports the results of the fits. It is clear 

2 

by looking at the that our data are not in 
the kinematical range where the chiral logarithms 
are visible. This is visualized in Fig. 2.a-b. A 
possible explanation for this behaviour is that, 
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08 - quenched logarithms 




08 - polynomial fit 



• Mj(M„,M.,E,)/FIT p=0 
«Mj(Mk,M,,E,)/FIT p=27i/L 
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Figure 2. (a)-(b) Quality of the qxPT and polynomial fits for Og. (c) xPT-polynomial matching for Og 
at = = 0.4 GeV, Mk = 0.41 GeV (only data with A/^ ^ Mk are plotted). 



since logarithms come from loops of Goldstone 
bosons (GBs), they become relevant at a scale 
of the order of the mass of the GBs. Since the 
contribution of the counterterms is instead con- 
trolled by the physics at the scale of the cut-off 
(i.e. 0(1 GeV)), it is dominant in the kinematical 
region accessible on the lattice (in our simulation 
M^,Mk,E^ e [0.54,1.2] GeV). 

On the other hand, one expects that predic- 
tions of xPT (including the logs.) are valid when 
all the energy scales are sufficiently small (let's 
say below the mass of the Kaon). If we assume 
that there exists a region in which the two de- 
scriptions of the form factors are reasonably close 
to each other, we can perform the matching be- 
tween these two forms [8]. This amounts to fit 
the data with (^(^•^**))p°'y and then to impose the 
equahty of {MgfY'"^ with (A1<If )''''^ and of 
their first derivatives with respect to the kine- 
matical variables M^, M^, p-^+-pK, Ptt+'PttO at a 
matching point (chosen in the SPQR kinematics). 

Results for different matching points are shown 
in Tab. 2. Lattice data have not yet been cor- 
rected for all of the finite volume corrections ex- 





mTJ2 




poly. 


0.714(62) 


0.120(11) 


match. 


XPT qxPT 


xPT qxPT 


0.3,0.5 
0.4,0.5 
0.3,0.31 
0.5,0.51 


0.697(60) 0.738(64) 
0.664(57) 0.758(66) 
0.638(55) 0.711(61) 
0.635(55) 0.825(72) 


0.117(11) 0.125(12) 
0.111(10) 0.128(12) 
0.106(10) 0.120(11) 
0.106(10) 0.141(13) 


Table 2 



Results (same units and scheme of Tab. 1) for the 
matching procedure at the point (A/jr = EjtjMk)- 



plained in [2] and so our results are still prelim- 
inary. Nevertheless it's interesting to study the 
systematics of the xPT-polynomial matching and 
the contribution of the NLO at the physical point. 
Since quenched logarithms are very different from 
the full QCD ones (and assuming that quenching 
does not change drastically the kinematical be- 
haviour of the form factors in the range accessible 
to the lattice) we think that the most sensible way 
of extrapolating to the physical point is to per- 
form the matching with full xPT. If we choose the 
matching point (Af^, )=(0.4,0.5)GeV we get: 

(Os) = 0.664(57) (1«) (Ifo) (Or) = 0.111(10) (t«) {H) 

where the first error is statistical, the second is 
the systematics of non-perturbative renormaliza- 
tion and the third is the systematics due to the 
variation of the matching point. At the physical 
point, the NLO gives a negative contribution of 
order 15 ^ 22% with respects to the LO. 
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